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RADIO PULSAR DEATH LINE REVISITED: IS PSR J2144-3933 ANOMALOUS? 
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ABSTRACT 

We reinvestigate the radio pulsar "death Unes" within the framework of two different types of polar 
cap acceleration models, i.e., the vacuum gap model and the space-charge- limited flow model, with either 
curvature radiation or inverse Compton scattering photons as the source of pairs. General relativistic 
frame-dragging is taken into account in both models. We find that the inverse Compton scattering 
induced space-charge-limited flow model can sustain strong pair production in some long-period pulsars, 
which allows the newly detected 8.5s pulsar PSR J2144-3933 to be radio loud, without assuming a special 
neutron star equation-of-state or ad hoc magnetic field configurations. 

Subject headings: pulsar: general - pulsar: individual: PSR J2144-3933 - radiation mechanism: 
non-thermal 



charge-limited flow model with frame-dragging included, 
but some degree of the point dipole offset is needed. 

2. DEATH LINES IN VARIOUS MODELS 

The death line has been defined (e.g. RS75) by the con- 
dition that the potential drop across the accelerator {AV) 
required to produce enough pairs per primary to screen 
out the parallel electric field is larger than the maximum 
potential drop ('I'max) available from the pulsar, in which 
case no secondary pairs would be produced 0. 

It is worth noting that pulsar death lines are actually 
model- dependent. Besides the mo del- dependent ^>max (see 
1^ and 0), the form of i^y within the accelerator (which 
depends on the boundary conditions and on whether the 
general relativistic frame-dragging effect is taken into ac- 
count), the typical energy of the 7-ray photons (which de- 
pends on whether their origin is curvature radiation (CR) 
or inverse Compton scattering (ICS)), and the strength 
of the perpendicular magnetic field the 7-ray photon en- 
counters (which depends on the field strength and the field 
line curvature near the neutron star surface) can change 
AV considerably and alter the death lines. Further- 
more, the obliquity of the pulsar (which changes and 
'I' max) and the equation-of-state of the neutron star (which 
changes the moment of inertia / = 10*^/45 and the radius 
R — 10^i?6 of the star) will also influence the location of 
the death lines. As a result, the phase space for pulsars to 
die should be a "valley" rather than a single line. CR93 
defined a death valley, within the framework of the RS75 
vacuum gap model, as the phase space range between the 
death line of the dipolar field configuration and the death 
line of some special multipolar field configurations. Here 
we will also adopt such a death valley, and regard it as also 
including the scatter of obliguities for different pulsars]^ 

^Laboratory of High Energy Astrophysics, NASA/Goddard Space Fhght Center, Greenbelt, MD 20771 
^National Research Council Research Associate Fellow 

^SM&A Corporation, Space Sciences Division, Upper Marlboro, MD 20774 

^Although some authors argued that pair production might not be an essential condition for pulsar radio emission (e.g. Weatherall & Eilek 
1997), no such pairless radio emission theory is fully constructed. 

^Note that there are some additional criteria to constrain the death lines for the millisecond pulsars (Rudak & Ritter 1994; Qiao & Zhang 
1996; Bjornsson 1996) and the pulsars with strong surface magnetic fields (Baring & Harding 1998), and we will not address them in this paper. 

^Strictly speaking, there could be different death lines for pulsars with different obliquities. The death lines could be quite different for the 
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1. INTRODUCTION 

Copious pair production in pulsar inner magnetospheres 
has long been conjectured as an essential condition for the 
radio emission of pulsars since the pioneering work of Stur- 
rock (1971)^ Generally speaking, pair production can cre- 
ate a dense plasma which may allow various coherent in- 
stabilities to grow, so as to account for the high brightness 
temperature observed from the pulsars. The secondary 
pairs (rather than primary particles) with 7 ^ 10^ — 10'* 
can produce typical radio-band emission within various 
models (e.g. Ruderman & Sutherland 1975, hereafter 
RS75; Melrose 1978; Qiao & Lin 1998). As a consequence, 
the so-called radio pulsar "death line" is defined as a line 
in a two-dimensional pulsar parameter phase space {P — P 
diagram, P — diagram, or P — $ diagram, where P is 
the pulsar period, P the pulsar spin-down rate, Bg the sur- 
face magnetic field, and $ the polar cap potential), which 
separates the pulsars which can support pair production 
in their inner magnetospheres from those which cannot 
(RS75; Arons & Scharlemann 1979, hereafter AS79; Chen 
& Ruderman 1993, hereafter CR93; Rudak & Ritter 1994; 
Qiao & Zhang 1996; Bjornsson 1996; WeatheraU & Eilek 
1997; Arons 1998, 2000). Present pulsar death fine theo- 
ries require some degree of anomalous field line configura- 
tions (multipole component or offset dipole) to interpret 
known pulsar data. However, a recently discovered long- 
period (8.5s) pulsar PSR J2144-3933 (Young, Manchester 
& Johnston 1999) is clearly located well beyond the con- 
ventional death valley (CR93), unless special neutron star 
equation-of-state or even ad hoc magnetic field configura- 
tions are assumed. This challenges the widely accepted be- 
lief that pair production is essential for pulsar radio emis- 
sion. Arons (1999) found that this pulsar is located within 
the death valley of an inverse Compton controlled space- 
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Modification of the equation-of-state or of the polar cap 
radius wiU also modify the death valleys systematically. 
The surface magnetic field for the star-centered dipolar 
configuration is Bd = 6.4 x 10i^(PP)i/2ji/2^^ 3 j-gg^rd- 
less of the internal field geometry (Shapiro & Teukolsky 
1983; Usov & Melrose 1995). Thus the only model depen- 
dence of the dipolar surface field is the offset of the dipole 
center from the star center. 

Two subgroups of inner gap models were proposed by 
adopting different boundary conditions. The vacuum- 
like gap model (V model, RS75) assumes strong bind- 
ing of ions at the neutron star surface {E^^{z = 0) / 0), 
while the space-charge-limited flow model (SCLF model, 
AS79) assumes free emission of particles from the surface 
= 0) = 0). Both of these models were originally 
proposed by assuming that CR of the primary particles 
is the main mechanism to create 7-ray seeds to ignite the 
pair production cascades, and by neglecting general rel- 
ativistic effects. They were improved later by different 
authors to include ICS and the inertial frame-dragging ef- 
fects. The role of frame-dragging in pulsar physics was 
first explored by Muslimov & Tsygan (1992, hereafter 
MT92) and updated by Muslimov & Harding (1997, here- 
after MH97) and Harding & Musfimov (1998, hereafter 
HM98) within the framework of the SCLF model. MT92 
noted that since stellar rotation actually drags the local 
inertial frame (LIE) to rotate with an angular velocity 
r^LiF — fl*Kg{R/r)^ (r2* is the angular velocity of the 
star), the electric field required to bring a charged particle 
into corotation is weaker than that in the flat spacetime, 
since such a fleld only needs to compensate the angular ve- 
locity difference between and JIlif rather than the dif- 
ference of J7* and the angular velocity at infinity (which is 
0). As a result, near the star surface, the Goldreich- Julian 
density become^ 
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0.27 (HM98; we 



will adopt a typical value of 0.15 hereafter), a = (1 



1/2 



0.78 is the redshift factor, is the gravita- 



tional radius, and M is the mass of the neutron star. 
2.1. Vacuum gap model (V model) 

In principle, pulsar i^n arises from the deviation of the 
local charge density (77) from the Goldreich- Julian density 
(77^). If the binding energy of the positive ions is large 
enough to prevent the ions from thermionic or field emis- 
sion ejection, a vacuum-like gap (RS75) will form right 
above the neutron star surface, in which the charge deple- 
tion is very large (of the order of 777? itself) so that a very 
strong .Ell is built up right above the surface. This picture 
was questioned since later calculations showed that the 
ion binding energy is actually not high enough (for recent 
reviews, see, e.g. Usov & Melrose 1995). Recent observa- 
tions indicate that some pulsars may favor the existence of 
such vacuum gaps (Vivekanand & Joshi 1999; Deshpande 
& Rankin 1999), and some ideas to solve the binding en- 
ergy problem have been proposed (e.g. Xu, Qiao & Zhang 

extreme cases of aligned and orthognal rotators. 

'^For an explicit expression for r;^, see eqs.[4], [8] and [16] of HM98. 



1999). The inclusion of ICS in such a model was carried 
out by Zhang & Qiao (1996) and Zhang et al. (1997), 
and the corresponding death line was examined by Qiao 
& Zhang (1996). The influence of the frame-dragging ef- 
fect on the model has not been examined before. 

The maximum potential available for the V model is just 
the homopolar generator potential, which is the potential 
difference between the pole and the edge of the polar cap 
at the star surface and reads 



*max(V) ~ {l-ng) 



2c2 



5.6xlO^^(Volts)Brf,i2P-^i?i 



2 733 
6- 

(2) 

Note frame-dragging modifies the flat spacetime result 
by a factor of C = (1 ^ '*g) ~ 0.85, which arises from 
(fj, — J^lif)- The solution of the one-dimensional Pois- 
son's equation (the infinitesimal gap in RS75) then gives 
Ei\{z) = C(2n^B/c){h- z) and AV = C(^^*B/c)/l^ which 
are essentially RS75's results with the correction factor (. 

The gap height depends on three length scales, i.e., 
/acc, the acceleration scale before the primary elec- 
tron or positron achieves a high enough energy 7c; 
le ^ c[y{'-fcmc^) / Ec]^^ , the mean free path of the elec- 
tron/positron with 7c to emit one 7-ray quanta with en- 
ergy Ec] and Iph = xpiBcri/Bs){2mc^ /Ec), the mean 
free path of the 7-photon before being absorbed, where 
Bcri = m?c' /he ~ 4.4 x 10^'^G is the critical mag- 
netic field, and x 0-1 is the key parameter to de- 
scribe 7 — i? absorption coefficient (Erber 1966). The 
electron mean free path 1^ should not exceed the gap 
height h = /acc + Iph- For the V model, one gets 
7c = 2({e/mc^){fl^B /c)la,ccih — Zacc/2) with the form 
of Eor the CR-induced cascade model, we have 

E, =. {3/2)ihc/phl and 1^ = i9/4){hc/e'')p-f-' « Iph - 
'acc ^ h. To treat the gap breakdown process, we are 
actually looking for the minumum of h, which could be 
obtained by setting the derivative h with respect to ^acc 
to zero. With some approximations, we finally get the 
gap parameters as h = 5.0 x 10^ {cm)(-^^'^ P^^"^ B'^^'^ pl^'^ , 

and AV = 1.6 x 10^^{Yo\ts)C^/'^ P-^/'^ B^^^'^ Pq^\ which 
are analogous to RS75 (their eqs.[22],[23] except for the ( 
correction), who treated the problem by simply adopting 
h ^ Iph- By setting AV — <i>max(V), we get the death fines 
of this model {( ^ 0.85 has been adopted) (cf. eqs.[6],[8] 
of CR93) 



logP= (11/4) log P- 14.62 
logP = (9/4) logP - 16.58 + logpe 



[I] 



(3) 

(4) 



for the dipolar and the p ^ R, Bs ^ Bd multipolar field 
configuration, respectively. 

For the resonant ICS-induced gap, we have Ec = 
2-fh{eB/mc) (Zhang et al. 1997; HM98), and jres ^ 
4.92 X 10"B^2767c"^ (Dermer 1990), so that 1^ ^ 
0.0027672p-^rg-\ Since Zacc < ~ Iph for the ICS 
case, we solve the gap height by setting h ^ Iph le- 
Also treating Tg self-consistently through self-sustained 
polar cap heating, i.e., T = {eAVN /aiTrp)^/^ , where 
N = [riB /ATre)TTrp is the polar cap luminosity, and Vp 
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is the polar cap radius, we finally get the gap height 
h = 2.6x 104(cm)pi/75-ii/7^4/7^_i/i4^ gj^p pQ^g^_ 

tial AV = 4.2 x 10^^(Volts)P-^/'^B;l^^'^ p^^'^C^/'^ , and the 
surface temperature Tg = 5.9Bi~2^^^P-3/7p^/7(-3/i4 (j^Q^g 
that such a high polar cap temperature conflicts with the 
observations). The death lines of this model are therj^ 

logP== (2/11) logP- 13.07 [II] (5) 
log P (-2/11) log P- 14.50 +(8/11) log p6 [II'] (6) 
for dipolar and multipolar configurations, respectively. 
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I, I': CR-V model 
il, If: ICS-V model 

III, 111'; CR-SCLF model 

IV, IV: ICS-SCLF model " 
PSR J3144-3933 

PSRs (Parkes MBS) _ 
Princeton PSR Catalog 



Log(P) (s) 



Fig. 1. — Death lines in various models. Solid lines are for dipole, 
and dashed lines are for p ~ lO^cm, Bs ~ -Bjj multipole. Other 
captions are marked in the plot. Pulsar data follow the latest results 
(Manchester et al. 1999). 



hues in this model are very high (see eq.[48], [50] of AS79). 
The ICS-induced version of such a SCLF model was pre- 
sented by Luo (1996). HM98 explicitly studied both CR- 
and ICS- controlled SCLF accelerators with the frame- 
dragging effect included, and such a model is very good 
in interpreting high energy radiation luminosities of the 
spin-powered pulsars (Zhang & Harding 2000, hereafter 
ZHOO). The inclusion of the general relativistic frame- 
dragging effect in such SCLF models (MT92; MH97) is 
essential in two respects. First, besides the B declination, 
rjji in curved spacetime has an extra (R/r)^ dependence 
(see second term in the bracket of eq.|l[), while the cur- 
rent conservation requirement leads to ry ~ — (1 — Kg) 
near the surface. As a result, E^^ is built up much faster. 
Secondly, the maximum potential available is much larger 
than the flat spacetime value (but smaller than that in V 
models), which is^ 



$^ax(SCLF) 



' 2c2 



1.0xl0^2(Volts)P<j, 
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2.2. Space- charge-limited flow model (SCLF model) 

If the charged particles (electrons or ions) can actually 
be pulled out freely from the neutron star surface (which 
is favored by ion binding energy calculations), a space- 
charge-limited flow is a natural picture, with Py = at 
the surface. The P|| at higher altitudes then arises from 
the small imbalance of the local charge density rj from rjn 
(eq.[||) due to the flow of the charged particles along field 
lines. The conservation of current requires cx r~^, thus 
the deviation (77 — rjn) arises from the extra r-dependence 
of Vr (besides B declination, which is oc r~^). In flat 
spacetime, this dependence is just the "flaring" of the 
field lines, on which the AS79's pioneering SCLF model 
is based. Since such a deviation is so small, it takes a long 
length scale for a particle to be accelerated to pair produc- 
ing energy via a gradual built-up of Py, so that the gap 
shape is usually narrow and long. The maximum poten- 
tial available in this model is much smaller than the one 
available in V models (sec footnote 9), so that the death 

*The typical ICS photon energy adopted here is the one for resonant 
the resonance with the photons at Planck's peak may become important 

^For a full expression of $inax(SCLF), see eq.[54] of HM98. What we 
term (the maximum available for the AS79 model). 



(7) 

As a result, the death lines are considerably lower than 
the AS79 model. 

By introducing an upper P|| = boundary at the pair 
formation front, HM98 have presented the explicit for- 
malism and detailed numerical simulation of P|| within 
SCLF accelerators. Unfortunately, simple analytic formu- 
lae applicable for all the cases are not available. How- 
ever, we notice that near the death lines, the height of 
the accelerators (h) are all larger than the polar cap ra- 
dius rpc = [n^R/cf{l)Y/^R (/(I) ~ 1.4 is the factor of 
curved spacetime), and that the Py in the accelerators has 
achieved saturation (see ZHOO for discussions of different 
regimes of acceleration Py and the criterion to separate 
the two regimes, their eq.[53]). In such a case, we can 
adopt the following approximate acceleration picture: near 
the star surface, Py grows approximately linearly with re- 
spect to the height z (thus is analogous to V model in this 
regime) with the form Py(l) — [3^^/(1 — e)]{^l^,B /c)z ~ 
l57Bi2P-^z (eq.A[3] of HM98), where e = rg/R ~ 0.4, 



and saturates above z ~ rpc/3 at a value of Py(2) 
{3Kg/2){n, ■B/c)rp,{rpjR){l~e) ^ 3.5 x lO^PiaP-^P^ 
(eq.[A5] of HM98), where /(I) = 1.4, ^ = 0.7 have been 
adopted. To study the death lines, we assume that the 
accelerators are located at the surface, though HM98 have 
argued that accelerators could be (0.5 — 1)P above the sur- 
face in young pulsars due to anisotropy of the upward ver- 
sus downward ICS. In old pulsars, the returning positron 
fraction is smaller, and the lower pair formation front may 
not exist, so that both the CR- and the ICS- induced SCLF 
accelerators could be formed at the surface. 

For the CR-induced model, we again have le «C lace ^ 
Iph- Since h > Vpc near the death lines (ZHOO), we can 
adopt 7c ~ (e/mc2)Py (2)Zacc- Following the same proce- 
dure as for the V model, but using AV = Py(2)/i (re- 
call the quadratic form of V model and note the differ- 
ence), we get /i = 3.3 X \^^{<z^)P^I'^B^:^ p)!"^ R^"^^"^ and 

scattering. In the cases of high temperatures, the scatterings above 
(Zhang et al. 1997), and the death lines could be lower. 

adopted here is the cos x term, which is much larger than the sin x 
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= 3.5x 10"(Volts)P~i/2pg/^i?6^^ Equating Al/ with 
*max(SCLF) (eq.0]), we get the death lines 

logP= (5/2) log P- 14.56 [III] (8) 

logP = 21ogP- 16.52 + logp6 [in'] (9) 

for dipolar and multipolar field configurations, respec- 
tively. 

For the resonant ICS-induced SCLF model, we again 
have lace <C ?e ~ Iph and h > r^c- This brings an important 
difference with the CR-induced case, that is, one should 

adopt the linear form, e.g., jc — {e/mc^) E\<^{l)dz 
to describe the acceleration phase, but adopt the saturated 
£^11 form to describe the final potential, i.e., ISV = E\\ (2)ft.. 
Since the SCLF model has a much lower charge deficit than 
the V model, the number of reversed positrons required to 
screen the field is only a factor of / ~ | (i9i?|| /dr)/{8'Kp^j ) | 
of the Goldreich- Julian density (AS79; ZHOO), so that one 
gets less polar cap heating, i.e., T = {eAV N f / aTrr^y^^ , 
in this model. For the saturated accelerators, this fac- 
tor is roughly / ~ 5.7 x 10~^P~"'^ near the surface 
(eq.[71] of ZHOO). A self consistent treatment of T fi- 
nally leads toh^ 9.7 x 104(cm)p4/i3 5-22/13^8/137^-2/13^ 

AV = 1.0 X 10"(Volts)p-22/i3B-2'/'V6^"i?6^^", and 
Te = O.llp-12/1351/13^2/13^6/13^ ^^^^ ^j^g ^g^^^ ^.^gg 

are 

logP = -(3/11) logP- 15.36 [IV](10) 

logP= -(7/11) logP- 16.79+ (8/11) log P6- [IV'](11) 

In this model, the polar cap temperature is sustained 
by the thermal energy released from the crust deposited 
by the reverse flow of positrons and high energy pho- 
tons. The primary beam fluctuations (typical timescale 
~ ft./c ~ 3 X IQ-^s) may result in discontinuous illumina- 
tion of the polar cap. However, because of the "inertia" 



of photon diffusion from the relatively deep layers to the 
photosphere (e.g. Eichler & Cheng 1989), this process is 
unlikely to affect the average polar cap temperature. For 
example, it takes 6 x 10-3r5-i'^s for the photons to 
diffuse upthrough the surface from a depth of lOOg cm^2^ 

3. CONCLUSION AND DISCUSSIONS 

The death lines of different models are plotted in Fig.l. 
A remarkable fact is, though PSR J2144-3933 is beyond 
the death valleys of both the CR- and ICS- induced V 
models and CR-induced SCLF model, it is well above the 
star-centered dipolar death line of the ICS-induced SCLF 
model. Thus one does not need to introduce a special 
neutron star equation-of-state or anomalous field config- 
urations at all to maintain strong pair formation in this 
pulsar. The ICS-SCLF death lines of Fig.l imply a very 
large phase space for radio emission in the long period 
regime, thought previously to be radio forbidden. In fact, 
Young et al. (1999) argued that the population of pulsars 
with parameters silimar to PSR J2144-3933 is very large, 
since the detectability of such pulsars is very small due to 
their small polar caps. We expect more pulsars to be de- 
tected in this region. A general trend in Fig.l is that death 
lines in the ICS-induced models have much flatter slopes 
than those in the CR-induced model. This arises from the 
very different P-, Bg- dependences of both Iph and Ig for 
the resonant ICS processes. The saturation of in long 
period pulsars is critical in lowering the death lines in the 
SCLF models relative to those of the V modclsj^ The 
difference is more prominant in ICS- induced models. For 
the ICS-SCLF model, our death lines are lower than the 
ones derived in Arons (2000). The difference might be due 
to the different polar cap temperature treatments. 

We thank the referee for good comments and sugges- 
tions and Demosthenes Kazanas, G.J.Qiao, R.X.Xu, and 
Z.Zheng for interesting discussions or helpful comments. 
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I'^V models also have By saturation, but only near the death lines when h ~ rpc (RS75). SCLF models, however, can have By saturation 
much farther away from the death lines and thus allows a greater distance for pair production before ^max is reached. 



